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Abstract 

We show that one- dimensional Bohmian mechanics is unique, in that, the Bohm trajectories are the only solutions 
that conserve total left (or right) probability. In Brandt et al., Phys. Lett. A, 249 (1998) 265-270, they define 
quantile motion — unique trajectories are solved by assuming that the total probability on each side of the particle 
is conserved. They argue that the quantile trajectories are identical to the Bohm trajectories. Their argument, 
however, fails to notice the gauge freedom in the definition of the quantum probability current. Our paper sidesteps 
this under-determinedness of the probability current. The one- dimensional probability conservation can be used for 
higher dimensional problems if the wave function is separable. Several examples are given using total left probability 
conservation, most notably, the two-slit experiment. 
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1. Introduction 

Bohm's formulation of quantum mechanics ad- 
mits particle trajectories [l '2 .'! . Given a solution 
tfj(x, t) = R(x, t)e lS ^ x -^/ h of Schrodinger's wave 
equation for real R(x,t) and S(x,t), Bohm defines 
the equation of motion for the particles as, 

mx = -W-V<3, (1) 

where the quantum potential is Q — —h 2 V 2 R/(2mR) 
In order that an ensemble of particle trajecto- 
ries satisfy the predictions of standard quantum 
mechanics, the particle velocities are further con- 
strained to obey the so called guidance law, 
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X = — VS. (2) 

m 

The guidance law guarantees the equivalence 
of the particle ensemble density p p and the quan- 
tum probability density p^, since both densities 
evolve according to a continuity equation. From 
Schrodinger's equation the quantum probability 
density has, 

and as a result of particle conservation, 

^ + V-(p p x) = 0. (4) 

Therefore, if p p (x,t = 0) = p^,(x,t = 0) then by 
the expressions above the two densities will remain 
identical. However, the definition of the guidance 
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law [and in turn the equation of motion, Eq. (1)] 
is not unique, since one can add a divergence-less 
gauge to the definition of the particle guidance law 
Eq. ((21), and the equivalence of the particle density 
p p and the probability density p^ is maintained as 
before. Hence, Bohm's formulation of quantum me- 
chanics is under-determined |4l5j . There has been 
some efforts to resolve this problem [6 7 8 9. Either 
the proposal doesn't work, or if it does, it must add 
additional physical assumptions to Bohm's theory. 
In Brandt et al. |10|llj . however, they give a non- 
physical reason for choosing Bohm's definition of 
the velocity field in one-dimension. They introduce 
the concept of quantile motion. A particle trajectory 
xp(t) is generated such that, 



cp(t) 



P 



p(x, t) dx = constant. 



(5) 



This is recognized also as the cumulative probability 
function (CPF) for the probability density p(x,t). 
Additionally, one can think of this expression as rep- 
resenting the total left probability during the par- 
ticle's trajectory. Note, since the total probability 
is conserved in time, one could have used the total 
right probability instead. 

Their argument for the equivalence of the quan- 
tile trajectories and Bohm's trajectories rests upon 
the similarity between the two expressions for the 
velocity fields, 



xp 



Jp 
P 



and 



X'Qohv, 



P ' 



(6) 



Where j p is the particle current, and is the quan- 
tum probability current. They assume that the two 
currents are equal, hence, the quantile trajectories 
and Bohm's trajectories will be identical (if the ini- 
tial conditions are the same) . The two currents, how- 
ever, can be related by a divergence-less gauge as 
shown above. Therefore, they fail to show that the 
quantile trajectories are in fact the Bohm trajecto- 
ries because of the gauge freedom in the definition 
of the two currents. 

Below we present a proof that sidesteps this 
under-determinedness of the currents by not in- 
troducing the current concept in the first place. 
Working directly from the quantile motion idea 
and Schrodinger's wave equation, we show that in- 
deed the unique quantile trajectories for quantum 
mechanics are the Bohm trajectories. This gives a 
non-physical reason for selecting Bohm's definition 



of the velocity field (in one-dimension, at least), 
namely, the conservation of the total left (or right) 
probability. 

Later, we show that the one-dimensional total left 
(or right) probability conservation can be used to 
generate Bohm trajectories for higher dimensional 
problems if the wave function is separable. Finally, 
we provide several comparisons of the quantile and 
Bohmian trajectories: harmonic oscillator, free par- 
ticle, the two-slit experiment, and a two-dimensional 
infinite square well for a separable wave function. 
In all cases, the quantile trajectories are identical to 
the Bohm trajectories. 

2. Unique Quantile Trajectories in Quantum 
Mechanics are Bohmian 



In Brandt et al. [10111] . they show that one can 
construct trajectories by requiring that the total left 
(or right) probability is conserved along a particular 
trajectory. In their paper, they work with the total 
right probability, 



Q 



p(x,t) dx. 



(7) 



x Q (t) 



But for the following we use the total left probability 
(noting that P + Q = 1 for all time) , 



X P (t) 



p 



p(x,t) dx. 



(8) 



This is also known as the cumulative probability 
function (CPF) for the probability density p(x,t). 
The CPF is one-to-one and monotonically increas- 
ing. At each time, there is only one x p that satisfies 
Eq. ijHJ) for a constant value of P. Therefore, there 
is a unique trajectory such that the total left prob- 
ability is conserved. Given that P is constant, 



dP 
~dt 



C P (t) 



dp 

~dt 



dx + p(xp, t)ip = 0. 



(9) 



Where we assumed that the density is zero at the 
lower boundary. At once this can be solved for the ve- 
locity field that yields unique trajectories that con- 
serve total left probability, 



1 



cp(t) 



xp 



p(x P ,t) 



dp 
dt dx - 



(10) 



2 



So far the discussion has been quite general and 
Eq. (|10[) is a definition for trajectories for any given 
density p(x, t), whether it be quantum or not. To ap- 
ply these quantile trajectories to quantum mechan- 
ics, we begin with Schrodinger's wave equation in 
one-dimension, 



dtp h 2 d 2 ^ 



(11) 



and the similar expression for the complex conjugate 
wave equation. Writing the complex wave functions 
in polar form, ip(x,t) = R(x,t)e lS ^ x ' t ^ h for real R 
and S, and following the typical derivation of the 
quantum continuity equation, we get that, 

P^TT- (I 2 ) 



dp dp dS 



dt dx mdx ' r mdx 2 
Notice, that we did not write the right hand side 
of this continuity equation as dj /dx where j = 
pdS/mdx — this would have introduced the undcr- 
determinedness of the probability current since 
there is no unique anti-derivative. Inserting Eq. (112| 
into Eq. (flO]). 



1 



X P {t) 



xp 



p{x P ,t) 



(dp dS d 2 S 



\ dx mdx ^ mdx 



dx, (13) 



and integrating by parts the second term of the in- 
tegrand (again assuming that the density is zero at 
the lower boundary), 

dS 



Xp 



mdx 



(14) 



Which is the one-dimensional Bohm velocity field 
of Eq. |2j). We conclude, therefore, that the unique 
one-dimensional quantile trajectories are in fact the 
Bohm trajectories in quantum mechanics. 

3. Extension to Higher Dimensions 

The extension of the quantile motion into higher 
dimensions is also discussed in Brandt et al. [10]. 
They show that instead of the total left (or right) 
probability being conserved in one-dimension, that 
in higher dimensions the probability is conserved 
inside a volume that is enclosed by a surface of 
Bohmian trajectories. This property, however, is not 
unique to Bohmian mechanics. Any velocity field x 
will conserve the probability inside the volume in 
configuration space since [12] , 



dp 
~db 



= — pV • x and 



dJ 
~dt 



= + JV • x, 



(15) 



where p is the probability density and J is the Jaco- 
bian such that the volume changes as dV(t) = JdVo. 
The probability inside this evolving volume is, 



Pi, 



pdV{t) = / P J dV . 



(16) 



Which implies that dP- ln /dt = 0. Therefore, any ve- 
locity field x will conserve the total probability in- 
side a volume that in enclosed by a surface of tra- 
jectories solved from x. This is in contrast to what 
was found in the one-dimensional case in $2] where 
there was a unique velocity field that satisfied the 
total left (or right) probability conservation. 

However, the quantile motion concept can be used 
to generate trajectories in higher dimensions if one 
uses the marginal distribution for each coordinate 
analogously to the cumulative probability distribu- 
tion (CPF) in one-dimension. Suppose the system 
can be described by N Cartesian coordinates, then 
the corresponding definition of Eq. ([5]) is, 



Xi{t) 

Pi = I Pi(xi,t) dxi 



for i = 1,2, 



,N, (17) 



where pi{xi, t) is the marginal distribution for the i- 
th coordinate. Again, we assume that the particles 
are conserved so, 



dp 
dt 



N a 

Y — 

^ dx, 

i—l 



(pii 



0. 



(18) 



Partially, integrating this continuity equation we 
have that, 



+ oo Xi(t) +o 



dp 

dt 



dxi . . . dx 



N 



-oo — oo — oo 

N +oo Xi(t) +oo 

3=1. 



-oo — oo — oo 



dXj 



dxi . . . dxN — 0. (19) 



Interchanging the partial derivative with respect to 
time and performing the ±oo integrations, and as- 
suming that the density is zero at ±cx), only the j = 
i integrals on the second term survive, 



*;(*) 



dpi 
-z— dxi 
dt 



-oo +OC 



J •• J P^i d x i • ■ ■ dx^i . . . dxN — 0. (20) 



-oo — oo 
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Then with the expression above, and differentiating 
Eq. (fl~7l) with respect to time, 



dt 



PiXi 



pii dx\ . . . dx^i . . . dxN-(21) 



Hence, the «-th coordinate, in general, doesn't con- 
serve total left probability of the marginal distribu- 
tion since ii could depend on the other coordinates. 
Suppose, however, that ii — ii(xi,t) (i.e. the mo- 
tion along the i-th coordinate is independent), then 
dPi/dt = 0, and the total left probability is con- 
served for the marginal distribution pi. 

In higher dimensional Bohmian problems the 
guidance law [Eq. ^j] for the i-th coordinate be- 
comes, 



x 



1 dS(xi,x 2 , ■ • . ,x N ;t) 



m dxi 
If the wave function is separable, then, 

tp = lpl(xi,t)lp 2 (x2,t) ■ ■■1p N (x N ,t). 



(22) 



(23) 



The probability density is also separable, p — 
Pi{x\, t)p2{x2, t) . . . pn{xn, t), and the phase be- 
comes S — Si(xi,t) + Sz{X2,t) + . . . + Sfq(x Nl t). 

Which implies that ii = (l/m)dSi(xi,t)/dxi, or 
that the velocity field for the i-th coordinate is 
independent of the other coordinates, for all i — 
1, 2, . . . , N. Hence, one can use the one-dimensional 
total left (or right) probability conservation de- 
scribed in independently for each coordinate, to 
generate higher-dimensional Bohmian trajectories 
for a separable wave function. 

4. Examples 

Below are several examples that compare the 
quantile trajectories to the Bohm trajectories. The 
first three examples are one-dimensional, while the 
last example is for a two-dimensional separable 
wave function. In each example, the wave func- 
tion is non-stationary so that dp/dt ^ 0. From the 
wave function, the probability density is computed, 
p(x,t) = \ip(x, t)\ 2 . The Bohmian trajectories can 
be numerically solved using Eq. (|14|) or more con- 
veniently from this alternative form for the velocity 
field 03], 



(24) 



h f ip* dip / dx — -ipdif)* / dx 
2mi \ p 




Fig. 1. (color online). Harmonic Oscillator with wave function 
as a superposition of the ground and first excited states. The 
quantile trajectory points (+) are shown superposed on the 
Bohmian trajectories. The graph is in naturalized units. 

The quantile trajectories are numerically com- 
puted by the method described in Appendix YK\ In 
the figures below, the quantile trajectory points (+) 
are shown against the Bohm trajectory. In all cases, 
quantile motion is able to reproduce the Bohm tra- 
jectories. 



4.1. Harmonic Oscillator 

In this example, the harmonic oscillator wave 
function was taken to be a superposition of the 
ground and first excited states, 




e e 



-iEot/H 



_ * 2 2x 
2a^/n a 



iEit/h 



(25) 



where a = y/H/muj, and Ej = hio{j + 1/2). Natu- 
ralized units were used so that h — 1, u> — 3, and 
m = 1. The range of the time was t — nAt 6 [0, 3] 
for n=l,2,3,..., and the size of each time step was 
At = 0.1. The cumulative probability function [the 
right hand side of Eq. (jS])] was approximated by a 
series of trapezoids (see Appendix [XJ each with a 
width of A:r = 0.2, and the position range was x S 
[—5, 5] (an area where the density was essentially 
non-zero). In Fig. ([1]), the quantile trajectory points 
(+) are plotted superposed on top of the correspond- 
ing Bohm trajectories. We see that the quantile tra- 
jectories are, in fact, the Bohm trajectories. Smaller 
trapezoid widths and time steps will produce more 
accurate and continuous trajectories. 



4 




Fig. 2. (color online). Free particle with an initial wave func- 
tion of a Gaussian centered around zero. The quantile trajec- 
tory points (+) are shown superposed on the Bohmian tra- 
jectories. Even with only six trajectories shown the spreading 
of the wave packet is evident. The plot is using naturalized 
units. 

4.2. Free Particle 



The wave function for the free particle (assumed 
to be Gaussian initially with width a) was taken to 
be, 



' 2 a \ p-ax 2 /[l+(2ihat/m)] 

#r,t)=( — 



yjl + (2ihat/m) 



(26) 



Naturalized units were used so that % = 1, m = 
1, and a = n/2. Again, t = nAt 6 [0,3] (n 
1 . 2,3,.. .) with time steps of At = 0.1. The trape- 
zoid widths (see Appendix [A"f were Ax = 0.2, while 
the range was x G [—5,5]. In Fig. @, the quantile 
trajectory points (+) are shown against the Bohm 
trajectories. Notice that the ensemble of trajectories 
depict the familiar spreading of the wave function. 

4.3. Two-Slit Experiment 

This two-slit example is from §5.1.2 in Holland [3]. 
At first this problem seems to be two-dimensional. 
However, the motion along the coordinate from the 
slits to the screen [x in Fig. ([5])] is assumed uniform, 
thus the probability density is one-dimensional and 
is only a function of y and t. To allow for easier calcu- 
lation we rescaled the experimental numbers given 



in Holland, so that h = 1, m = 1, and t„ 



100 



(the time between the slits and the screen), and then 
rescaled the results back to the actual numbers. Us- 
ing the trapezoid method as described in Appendix 
the time t = nAt e [0, t max ] (n = 1,2,3,.. .) with 
a time step of At = 2.5 (1/40-th the total time). 
At each time, the cumulative probability function 



(this time a function of y) was approximated by a 
series of trapezoids of width Ay — 3.24169 (1/80- 
th the range of y). The range of y was restricted 
to a width between ±129.668 where the probabil- 
ity density was essentially non-zero. In Fig. jSJ), the 
quantile trajectory points (+) are plotted along with 
the Bohm trajectories. The ensemble of trajectories 
makes the familiar two-slit intensity pattern on the 
screen (located on the right hand side of the figure). 
The quantile trajectories match the Bohm trajecto- 
ries quite well, even in those regions where the prob- 
ability density is very close to zero (between the high 
intensity bands). 

4.4. Separable Wave Function in a 
Two- Dimensional Infinite Square Well 

In Fig. ([4]) is a comparison of the quantile 
trajectories and their Bohm counterpart for the 
two-dimensional infinite square well. The separa- 
ble wave function was assumed to be ip(x,y,t) = 
xp x (x,t)ip v (y,t), where, 



4>x 



•(«»*) = Vi ( sin (f ) 



2ttx 
sm | —j— | i 



-HExt/h 



(27) 



and a similar expression for ip y (y,t). The energy 
Ei = tt 2 H 2 /2mL 2 , and naturalized units were used 
so that m = 1, H = 1, and the width of the well in 
each direction taken to be L = 1. The time was t — 
nAt £ [0, 1] for n = 1,2,3..., and the size of each 
time step was At = 0.05. The (x(t),y(t)) position 
of each particle was computed by approximating 
the cumulative probability function for each coordi- 
nate's marginal distribution by a series of trapezoids 
(see Appendix lA]) of width Ax = Ay = L/30. In 
Fig. (J4j) , the quantile trajectory points (+) arc plot- 
ted superposed on top of the corresponding Bohm 
trajectories. For the separable wave function, the 
quantile trajectories are again identical to the Bohm 
trajectories. 

5. Concluding Remarks 

We showed that for one-dimensional probability 
densities one can define unique trajectories that con- 
serve total left (or right) probability. These trajec- 
tories for one-dimensional quantum mechanics are 
the Bohm trajectories. What is interesting is that 
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Fig. 3. (color online). Two-Slit Experiment as described in §5.1.2 of P.R. Holland, The Quantum Theory of Motion: An Account 
of the de Broglie-Bohm Causal Interpretation of Quantum Mechanics, (Cambridge University Press, New York, 1993). The 
quantile trajectory points (+) are shown superposed on the Bohmian trajectories. The initial positions in each slit (left side 
of figure) are assumed to be Gaussian, and the ensemble of trajectories makes the familiar bands of bright and dark on the 
screen (right side of the figure). 




0.0 — 
0.0 



0.4 



0.6 



0.S 



one can generate these one-dimensional Bohm tra- 
jectories without introducing any additional physi- 
cal concepts (for example, the quantum potential). 
Hence, most of the example calculations showing 
Bohm trajectories that have appeared in the liter- 
ature over the years, can be generated without any 
additional physical assumptions. In addition, using 
the quantile motion idea one can generate Bohm-like 
trajectories for any density, whether or not a quan- 
tum wave function is known or is exists. In fact, the 
density could have been found experimentally, and 
the trajectories computed that depict its the evolu- 
tion. The quantile motion idea can be used to gen- 
erate higher dimensional trajectories for separable 
wave functions. 



x(t) 



Fig. 4. (color online). Comparison of the quantile trajectories 
(+) and the Bohm trajectories for a separable wave function 
in the two-dimensional infinite square well of size 1 X 1 in 
naturalized units. The initial position of each trajectory lies 
on the line from (0.5,0) to (0,0.5). 
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CPF(x) 




Fig. A.l. At each time step the cumulative probability func- 
tion (CPF) curve (in grey) is approximated by a series of 
trapezoids of width Ax. The position xp, corresponding to 
the constant quantile P value, is found by solving the linear 
equation of the top line of the particular trapezoid. 



Appendix A. Trapezoid Method for 
Generating Bohmian Trajectories Prom 
Probability Conservation 

The method described below to generate one- 
dimensional Bohmian trajectories does not use 
either of Bohm's equations of motion, Eqs. |U [5]). 
Instead, the trajectories, xp(t), are computed using 
the quantile motion concept (or the conservation 
of the total left (right) probability) in Eq. for 
particular constant P values. This expression, in 
general, can not be solved in closed form, and must 
be solved numerically. 

The trapezoid method approximates the cumu- 
lative probability function (CPF) with a series of 
equal width, Ax, trapezoids, see Fig. (|A.ljl . At each 
time t = nAt (n = 1,2,3...) with time step At, the 
height of each side of every trapezoid is found by nu- 
merically integrating Eq. ([8]). A particular quantile 
trajectory has a constant P value between zero and 
one. From the P value for a particular trajectory 
the corresponding trapezoid is found. The xp value 
at each time is then calculated by solving the linear 
equation for the top segment of the corresponding 
trapezoid. As the number of trapezoids used for each 
time step increases, the better the approximation of 
the CPF curve, and hence the better the calculated 
position of the particle at the specific time step. And 
as the size of the time step is reduced the smoother 
the resulting trajectories. We found for the exam- 
ples in $4] that Ax/ At ~ 3 gave fair results. 
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